Generation of lattice structures of optical vortices 
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We demonstrate experimentally the generation of square and hexagonal lattices of optical vortices 
and reveal their propagation in a saturable nonlinear medium. If the topological charges of the 
vortices are of the same sign the lattice exhibit rotation, while if alternative, we observe stable 
propagation of the structures. In the nonlinear medium the lattices induce periodic modulation of 
the refractive index. Diffraction of a probe beam by this nonlinearity-induced periodic structure is 
observed. 



I. INTRODUCTION 

Optical vortices are intriguing objects that attract 
much attention (!]] and display fascinating properties 
with possible applications in the optical transmission of 
information, or guiding and trapping of particles. They 
have a characteristic screw-type phase dislocation [||, 
which order multiplied by its sign is referred to as a topo- 
logical charge (TC). The study of optical vortices and 
more general - phase singularities, suggests not only new 
directions of fundamental research but also provides links 
to other branches of physics, as quantum optics su- 
perfluidity [Q, Bose-Einstein condensates j^H, and cos- 
mology. 

Optical vortices can be generated in several different 
controllable ways: in lasers with large Fresnel numbers 
[7|, by helical phase plates |^|, laser mode converters 
[9,|l0[, or computer-generated holograms (CGHs) JlT[ . 
The method of the CGHs, however, is the most com- 
monly used, since it allows precise control of the vortex 
position, TC, and possibility for generation of specific 
patterns of optical vortices. 

The propagation dynamics of a single vortex, both in 
linear and nonlinear media, has been a subject of many 
researches (see, e.g., [|L2|-|l4j]), whereby also the non- 
canonical properties of the vortex have been taken into 
account jl5],[l6| . It has been shown that the vortex posi- 
tion on the background beam is strongly affected by any 
source of phase and intensity gradients fll7j-|l9[ and could 
be controlled by interference with a weak plane- wave p0| . 
Special attention attract the vortices propagating in a 



self-defocusing nonlinear medium (NLM) , where they can 
form an optical vortex soliton (OVS) pi] . (For overview 
on OVSs see Ref. |j], chapters 7-8.) OVSs induce in 
the medium optical waveguides 1 22-24] which can guide 
weak information beams. An OVS was first experimen- 
tally generated in Kerr NLM |2[| and later in media with 
other types of nonlinearity: saturable-atomic pq , pho- 
torefractive ^] , and photovoltaic Q . Recently an OVS 
was also observed in a quadratic NLM with defocusing 
response. However, attention has been paid to avoid the 
modulational instability of the plane-wave background 
beam |§. 

The propagation of multiple-charged OVSs has also 
been investigated |^,^lj. It was found that they are 
topologically unstable and decay into vortices of unit 
charge |32j . The vortices produced by the decay, can ar- 
range themselves in regular patterns (vortex ensembles) 
while interacting with each other by phase and intensity 
gradients. The decay of the higher-order vortices obey 
the general principal of conservation the total angular 
momentum (AM) of the beam carrying them. Addition- 
ally, for a closed region of space the net topological charge 
must be conserved under continuous evolution provided 
that no vortices enter or leave the region. 

An ensemble of optical vortices exhibits a fluid-like mo- 
tion (l^]3^] which strongly depends on the geometrical 
configuration. The propagation of the simplest vortex 
ensemble, namely a vortex pair, has been investigated by 
several groups fh|[l|||||5) . In Ref. Q the rotation of 
the pair of vortices with equal TCs is controlled by the 
Gyou phase of the host Gaussian beam. By changing 
the beam intensity, the position of the beam waist inside 
the self-defocusing NLM changes, thus changing the an- 
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gle of the rotation at the output plane. A comparison 
between the degree of rotation of a vortex pair in linear 
and nonlinear regime was also performed in Ref. f35|. It 
was pointed out that the effect of rotation in the nonlin- 
ear regime could be more than three times higher than 
in the linear one. The enhancement is assigned to the 
nonlinear confinement of the vortex cores, which allows 
the vortices to propagate as vortex filaments. 

Recently the propagation of vortex arrays has been in- 
vestigated. Such arrays were generated by a bent glass 
plate j36|, or as a result of transverse instability of dark 
soliton stripes [p7|-p9[ . The instability could be enhanced 
additionally when the dark-soliton stripe interacts with 
an optical vortex, causing "unzipping" of the stripe jh^] . 
Ensembles of ordered optical vortices have also been in- 
vestigated in quadratic NLM and are promising for con- 
trollable generation of multiple- vortex patterns (4^] . The 
proposed method paves a way for creation of recon- 
figurable vortex ensembles by means of seeded second- 
harmonic generation. 

Regarding the fluid-like motion of the vortex- 
ensembles, a stationary configuration of vortices was 
found HI . It consists of three vortices of equal TC sit- 
uated in an equilateral triangle and an additional vortex 
with opposite TC in the center. That configuration was 
proven to be stable under small displacement of one of 
the dislocations. However, if the vortices are of higher- 
order they decay and subsequently form another station- 
ary configuration resembling a part of hexagonal honey- 
comb lattice. This fact directs attention to the investiga- 
tion of optical vortex lattices and characterisation of the 
propagation of the beams they are imposed on. 

Up to now lattices of optical vortices propagating in 
NLM were considered only theoretically. The simplest 
case of square lattice consisting of vortices with alter- 
native charges was investigated by direct modeling of 
four of them under periodic boundary conditions M2| . 
Later, lattices with different geometries imposed on a fi- 
nite background beam (conditions closer to experimental 
ones) were considered Q. It was shown that depend- 
ing on the TCs the vortex lattices can exhibit rotation 
or rigid propagation for equal or alternative TCs, respec- 
tively. In addition, lattices possess elasticity against dis- 
placement of one or more vortices out of their equilibrium 
positions. 

Here we report the first to our knowledge experimen- 
tal investigation of lattice structures of optical vortices 
in self-defocusing NLM. We concentrate our attention 
on two types of lattice geometries - square and hexag- 
onal one. When propagating in the NLM they induce 
a periodic modulation of its refractive index. For high 
beam intensities these changes are sufficient to cause a 
diffraction of a probe beam propagating perpendicularly 
to the volume with periodically modulated refractive in- 
dex. This diffraction may be controlled by steering the 
propagation of the vortex-lattice, e.g. by controlling its 



degree of rotation (for lattices consisting of vortices with 
equal charges). Additional control may be attained by 
changing the pump beam intensity which changes the re- 
fractive index of the medium and therefore the diffraction 
efficiency of the induced periodic phase grating. 

The maximal refractive index change in our experi- 
ment is of the order of 10~ 4 to 10~ 3 which is not enough 
to form an effective two-dimensional photonic band-gap 
structure p3[ . As a proof-of-principle, however, one can 
consider the possibility to trap glass spheres Q] by the 
optical vortices ordered in a lattice. This might give an 
opportunity for generation of effective two-dimensional 
photonic crystals. Such a crystal could be reconfigured 
by altering the degree of rotation of the lattice (by chang- 
ing the intensity of the focused background beam as in 
Ref. |54|) for equal TCs, or by use of dynamically recon- 
figurable holograms [ ff5f . 

We would also like to emphasize the close link between 
our results and the field of Bose-Einstcin condensates, 
where experimental investigations on vortex ensembles 
p6| , vortex arrays as a result of dark soliton-stripe in- 
stability 1 47 1, and vortex lattices have been reported 
recently. 



II. GENERAL ANALYSIS 

Let us consider the propagation of a beam in self- 
defocusing NLM with saturable nonlincarity where its 
evolution is described by the normalized nonlinear 
Schrodinger equation (NLSE) for the slowly-varying am- 
plitude envelope 



dE 1 , 



\E\ 2 E 



{l + s\E\2p 



0, 



(1) 



were Aj_ is the transverse Laplace operator. The trans- 
verse coordinates (x, y) are normalized to the character- 
istic size of the dark structures a, and the propagation 
coordinate z is normalized to the diffraction length of the 
dark beams. The background beam intensity / = \E\ 2 
is expressed in units of the intensity necessary to form 
one-dimensional (ID) dark soliton I\Dsoi of size a. The 
saturation parameter is defined by s = Imsoi/ 1 sat, where 
I sa t is the saturation intensity retrieved by the experi- 
mental conditions. The model of the saturation we use 
is introduced phenomenologically in order to describe 
the nonlinear response of the thermal medium. It was 
derived from a test experiment for self-bending of the 
background beam and described in detail in our previous 
works ||^,|^] . The parameters of the nonlinear response 
function s and 7 depend on the particular realization of 
the experiment, e.g. the properties of the NLM and the 
focusing conditions. In all measurements reported here 
we use thermal nonlinearity, and in particular ethylene- 
glycol dyed with DODCI (Diethyloxadicarbocyanine io- 
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dide). Two concentrations of the dye were used so that 
s = 0.4 and 1.2, respectively and 7 ~ 3 in both cases. 

In order to investigate the propagation dynamics of 
vortex lattices we first conducted numerical simulations 
by use of beam propagation method. The initial condi- 
tions were modelled as superposition of vortices situated 
in the nodes of a lattice: 
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with a square and hexagonal symmetry respectively. In 
Eq. (||) fjk are the nodes of the Bravais lattice repre- 
senting the physical lattice structure. The square lattice 
(Fig. 1 - upper row) coincides with its Bravais lattice, 
however the hexagonal honey-comb lattice (Fig. 1 - sec- 
ond row) is represented by as Bravais lattice with a base 
containing two vortices. If one define the primitive vec- 
tors of the Bravais lattice as b and c then its nodes are 
described as fjk — jb + /cc, with j, k integer numbers. 
The primitive vectors of the square lattice are orthogo- 
nal to each other and can be expressed in (x, y) coordi- 
nates as b = (A, 0), c— (0, A), where A is the distance 
between two neighbouring vortices. For the honey-comb 
lattice the primitive vectors are not orthogonal and are 
(\/3A,0) and c= (#A, |A). Then 



expressed as: b = (V3A,UJ and c 
the two vortices inside the elementary cell have positions 
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FIG. 1. Background beams containing vortex lattices of 
different geometries, (a) at the input of the NLM and (b,c) 
at z = 10 for lattice with equal and alternating TCs, re- 
spectively. Upper row: Images of square-shaped lattice; Sec- 
ond row - hexagonal lattice; Bottom row - transverse slices 
of hexagonal lattices. For comparison, transverse slice of the 
background beam (without vortices nested in) is shown with 
dashed curve. In all cases A = 5.0, 7o = 1, and s = 0.4. 
The images are grayscale-coded and white corresponds to the 
maximal intensity. 



The functions sq(r — fjk) and hex{f—fjk) describe the 
structure of the elementary cell of the Bravais lattice and 
are expressed as: 



sq(r - fjk) = tanh(|f - fjk\) exp 
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The sign function sgn is equal to +1 for equal TCs and 
— 1 for alternative ones. 

Then the lattice structure is imposed on a super- 
Gaussian (flat-top) background beam 



B(x, y, z = 0) = V A) exp < - 
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where the width w is chosen to exceed the characteristic 
width of the dark structures a more than 40 times, and 
Io is the maximal background beam intensity. 

We modelled the propagation of lattices of different ge- 
ometries and different TC distributions (see Fig. 1). No 
qualitative differences were observed in the propagation 
of vortex-structures with respect to the lattice geometry 
(square or hexagonal). The mode of propagation, how- 
ever, crucially depends on the vortex charge distribution 
(equal - Fig. lb or alternative - Fig. lc). Two character- 
istic differences are clearly seen: (i) In the case of equal 
TCs (sgn = +1) the superposition of the phases of all 
vortices results in azimuthal phase gradient and nonzero 
total AM which causes rotation of the whole structure 
(Fig. lb). In the case of alternative TCs (sgn = —1) 
the superposition of all the phases gives, in average, no 
phase gradient and zero total AM. As a result, steady 
propagation of the lattice is observed in the simulations 
(Fig. lc); (ii) In the case of equal TCs the non-zero total 
AM and the centrifugal forces lead to increased broad- 
ening of the background beam. The maximal intensity 
rapidly decreases along the NLM (J ~ 0.6 at z = 10 
Fig. lb - bottom row). The dependence on the intensity 
of the background beam in this case is relatively weak 
and the topological effects dominate the nonlinear ones. 
In the case of alternative TCs (Fig. lc) the background 
beam broadening is an effect only due to the combined 
action of the diffraction and self-defocusing nonlinearity 
and depends on the intensity. 

The degree of rotation of both lattice geometries for 
the case of equal TCs is depicted in Fig. 2. The rota- 
tion is due to the phase gradient, which is higher for the 
denser (square-shaped) structure. Therefore its rotation 
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(open squares in Fig. 2) is faster than the rotation of the 
hexagonal one (solid circles) . The dependence on the dis- 
tance is not linear because in the course of propagation 
the background beam spreads out and the distance be- 
tween the vortices increases. That causes a decrease of 
the angular velocity with the propagation length. 



neutral filters (F) are used to avoid its saturation. 




4 6 S 

Propagation distance, z 
FIG. 2. Angle of rotation of vortex lattices of equal TCs vs. 
propagation distance for square (open squares) and hexagonal 
(solid circles) geometry. The solid lines are drown to guide 
the eye. The lattice parameters are the same as in Fig. 1. 

Similar behaviour was previously described for the case 
of Kerr nonlinearity |3j| . Here we would like to point out 
the features caused by the saturation of the nonlinearity. 
As already mentioned, the rotation of the lattice and the 
increased beam spreading in the case of equal TCs are 
topological effects. The effects which depend on the non- 
linearity are related to the local intensity via the specific 
beam shape and the vortex-pattern formed on the back- 
ground. For example, the transverse profile of a single 
OVS in saturable medium differs significantly at differ- 
ent values of the saturation parameter |2(| (the OVS is 
broader at higher saturation). In the case of periodi- 
cally ordered vortices when each individual dark beam 
starts to broaden, the overlapping with the wings of its 
neighbours increases, while the individual cores do not 
change significantly. Since the vortices are imposed on 
a finite background beam and its total energy is con- 
served, bright peaks form in between as a result of local 
energy redistribution (See Fig. Ic - bottom row) . There- 
fore, even in saturable medium the well defined periodic 
modulation of the refractive index of the medium is still 
preserved. 



III. EXPERIMENTAL INVESTIGATION 

The experimental setup is similar to the one used in 
our previous works ]3^ , ^9| ] and is shown in Fig. 3. The 
488nm line from an Ar + laser is used to reconstruct the 
photolithographically produced CGH with the desired 
vortex-lattice. The +1 (or —1) order of the diffraction 
is separated by an iris diaphragm (D) and is focused on 
the input face of a glass cell containing the NLM. The 
output face of the cell is imaged to a CCD camera, and 
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FIG. 3. Experimental setup 
hologram; D - iris diaphragm; 
7.0cm and 8.0cm, respectively; 
- nonlinear medium; F - neutral filters; CCD 
characteristic distances are shown. 



| CCD | 

CGH - computer generated 
Li, L,2 lenses of focal length 
M - mirror; E - screen; NLM 
camera. The 



To assure the correct generation of the lattices by the 
CGHs we first opened the diaphragm allowing the +1 
diffraction order to interfere with the plane zeroth one. 
Interference patterns for three different cases are pre- 
sented in Fig. 4. The vortices appear as forks of in- 
terference lines. Two neighbouring vortices at each im- 
age are marked with arrows. The images show correctly 
reproduced square-shaped lattice with alternative TCs 
and two hexagonal lattices with alternative and equal 
TCs, respectively (Fig. 4, left to right). The images are 
brighter in the right-hand side since they overlap inho- 
mogeneously with the zeroth order beam. That inho- 
mogeneity also introduces an intensity gradient in the 
structure of vortices, which causes shrinking and vortex 
displacement from their regular positions. In overall that 
leads to deformation of the lattice. Being aware of this 
fact in the experiment we preserved the regular lattice 
structure by placing the diaphragm as close as possible 
to the CGH. 

square lattice hexagonal lattices 




alternative TCs alternative TCs equal TCs 

FIG. 4. Interferograms of three experimentally generated 
lattices. From left to right: square lattice with alternative 
TCs; hexagonal lattice with alternative TCs; and hexagonal 
lattice with equal TCs. Two neighbouring vortices in each 
image are marked with arrows. 

The features of the nonlinear propagation are deter- 
mined by measuring the characteristic nonlinear param- 
eters of the media for the two concentrations of the dye. 
For the lower one, the power necessary to form a ID 
dark-soliton stripe was estimated to be P\Dsoi — 22mW 
and the saturation power was P sa t — 60mW (measured 
in a self-bending scheme) [^9|. For the higher concen- 
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tration the characteristic powers were P\Dsoi — 20mW 
and P sa t — 16mW. The intensity distributions for two 
hexagonal lattices (with alternative and equal charges) 
at the end of the NLM with lower dye concentration are 
shown in Fig. 5. Because of some technical restrictions 
in synthesizing the CGH, for the lattice with equal TCs, 
the number of vortices encoded is less than in the holo- 
gram with alternative TCs. The geometrical characteris- 
tics (the elementary cell of the lattices), however, are the 
same in both cases. The propagation behaviour for both 
lattices is clearly different. While the lattice with alter- 
native charges exhibits steady propagation (Fig. 5a) the 
one with equal ones (Fig. 5b) tends to rotate (at about 
28° counter-clockwise). The background beam spreads 
stronger than in the case of a lattice with alternative 
charges. Unfortunately due to the different number of 
vortices, this fact is not obvious in Fig. 5. The smaller 
number of vortices in Fig. 5b modulates the background 
beam in a way that more filters were used to avoid the 
saturation of the CCD camera. As a consequence, the 
wings of the background beam in Fig. 5b are not seen 
and the beam diameter seems to be smaller as compared 
to Fig. 5a. In order to illustrate that the spreading is 
indeed higher in the case of equal TCs we looked in de- 
tail on the size of the elementary cell of the honey-comb 
lattice. Since the distances between the neighbouring 
vortices were encoded in the CGHs to be the same (the 
produced holograms were inspected by a microscope) any 
difference in the vortex separation is due to the evolution 
during propagation. In each image in Fig. 5 we inset 
the exact size and orientation of the elementary hexag- 
onal cell of the lattice (see bottom-right corner of each 
image). Indeed the comparison between the elementary 
cells of the lattices in both cases shows 18% bigger size 
for the one with equal TCs. 
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FIG. 5. Experimental images of the vortex lattices after 
10cm propagation in the NLM. (a) - hexagonal lattice with 
alternating TCs for powers of 10m W and 50mW. (b) - hexag- 
onal lattice with equal TCs for the same powers. The insets 
in each image represent the size and the orientation of the 
elementary cell of each lattice. 



The influence of the nonlinearity can be seen if the 
corresponding images for two particular powers are com- 
pared. We note again that the increase of the beam power 
does not influence the degree of rotation of the lattice 
presented in Fig. 5b since the waist of the laser beam 
is near the input face of the NLM. The higher power of 
the laser beam contributes, however, to the background 
beam broadening. Comparing the size of the elemen- 
tary cell of the same lattice at two different powers we 
estimated 15% broadening of the beam for the case in 
Fig. 5a and 12% for the case in Fig. 5b. That difference 
we attribute to the increased background beam size at 
the entrance of the NLM, which is due to the topological 
interaction of the equally charged vortices between the 
CGH and the NLM. 

The square-shaped lattices were investigated in the 
same way and qualitatively similar features were ob- 
served. We also investigated lattices with intentionally 
encoded defects in their structure, e.g. when one of the 
vortices is missing or all of the vortices in a line are shifted 
out of their equilibrium position. These experiments re- 
vealed the interesting property that the lattices exhibit 
elasticity. However the resolution in our experiments was 
not sufficient to resolve this feature in more than quali- 
tative manner. 

IV. DIFFRACTION OF A PROBE BEAM BY 
VORTEX LATTICES 

When an intense laser beam propagates along a NLM 
its refractive index changes proportionally to the inten- 
sity distribution. Since the vortex-lattices posses a peri- 
odic intensity distribution (see Fig. 1 and Fig. 5) one can 
expect periodic modulation of the refractive index. In 
a self-defocusing medium the higher local intensity will 
correspond to lower local refractive index. The lattices 
are imposed on a finite background beam and as whole 
it induces in the NLM a cylindrical defocusing lens (con- 
sidered in a perpendicular direction) , which is modulated 
by the dark vortex structure composing the lattice. In 
a thermal NLM, as slightly absorbing liquid, the nonlo- 
cal effect coming from the heat transfer also influence the 
refractive index change and effectively decreases its mod- 
ulation. The nonlocality is not taken into account in the 
model of Eq. (0). Its main influence is that at zero inten- 
sity (the points of vortex-phase dislocation) the refractive 
index change is non-zero (see a description in Ref. p(i|]). 

To investigate the modulation of the refractive index 
in the NLM caused by the presence of lattices we con- 
ducted an experiment in which a (probe) single-mode 
He-Ne laser beam was directed perpendicularly to then 
(pump) Ar + laser beam, as shown in Fig. 3. We aligned 
the probe beam in a way to cross the pump one 1cm in- 
side the NLM. (The higher dye concentration was used 
in this experiment). The input profile of the He-Ne laser 
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beam is shown in Fig. 6a and its circular symmetry is 
evident. When it crosses the Ar + laser beam the sym- 
metry is distorted and the beam elongates in direction 
perpendicular to the plane of Fig. 3. 




-4 



-5 



6 HIBH 

cm ^■I^H 

a) b) c) d) e) 

FIG. 6. Images of the probe He-Ne beam on the screen E. 
(a) the input He-Ne laser beam; (b) the intensity profile of the 
probe beam at low pump power (lOmW) for different parallel 
vertical displacements with respect to the pump one; (c), (d) 
diffraction of the He-Ne beam from the periodic phase grat- 
ings induced in the NLM by square and hexagonal lattices, 
respectively (pump power 80m W); (e) the diffraction pattern 
when a single vortex is imposed on the pump beam (pump 
power 30m W). 

First we identified the effect of optically-induced Gaus- 
sian cylindrical lens on the probe beam. In order to keep 
the same background beam characteristics the CGH was 
shifted in a way that only a region with parallel inter- 
ference lines was illuminated thus ensuring unmodulated 
background beam. This unmodulated pump beam in- 
duces a cylindrical lens in the NLM whereas the probe 
one passes trough the lens and defocuses (Fig. 6b - central 
image) . The diameters of both beams were estimated to 
be approximately equal at the cross point. Therefore one 



should expect that the probe beam will be strongly af- 
fected by aberrations of the induced lens. In Fig. 6b five 
probe beam profiles are shown for different positions of 
the He-Ne laser beam. The pump power is kept lOmW. 
Different input positions of the probe beam are achieved 
by parallel vertical translation by a simple periscopic sys- 
tem denoted for simplicity as mirror M in Fig. 3. The He- 
Ne beam elongates symmetrically if it crosses the pump 
in the center and asymmetrically if it is shifted up or 
down. At higher powers the aberration of the induced 
cylindrical lens becomes vertically asymmetric, probably 
due to the asymmetric heat diffusion in the cell. 

The situation is different when the vortex-lattice is im- 
posed on the background beam. At power of the Ar + 
laser beam higher than 20mW the vortices have well con- 
fined cores. Due to the nonlinear change of the refrac- 
tive index the vortex lattice 'writes' a phase grating in 
the NLM. The perpendicularly propagating He-Ne laser 
beam passes trough this grating and develops well pro- 
nounced diffraction orders at the output screen as seen 
in Fig. 6c, d at 80mW pump power. The constant of 
the phase grating written is apparently different for the 
square-shaped lattice (Fig. 6c) and for the hexagonal- 
shaped one (Fig. 6d). In the first case the period of the 
vortex-structure is smaller (so this of the phase grating) 
and higher dispersion in the diffraction orders is observed 
(higher angle of diffraction) . At lower powers diffraction 
orders were also observed. However, it was more difficult 
to distinguish them at the screen (E) because the effec- 
tive cylindrical lens had larger focal length. At different 
powers the magnitude of the refractive index changes and 
the modulation depth of the phase grating written in the 
NLM is different. This influences the energy redistribu- 
tion between the diffraction orders. Moreover, because 
of the finite number of the vortices in the lattices and 
the nearly equal sizes of the pump and the probe beams 
the diffraction from the phase grating couldn't be com- 
pared directly with the diffraction from an infinite pe- 
riodic structure. In our opinion the ratio between the 
intensities of the different diffraction orders is gradually 
influenced by the fact that different parts of the probe 
beam pass trough different number of vortices. Further, 
at the exit of the phase grating, the modulated probe 
beam is additionally affected by the aberration of the 
thermal lens. 

To assure that the observed diffraction structure is re- 
ally induced by the periodicity of the vortex lattices we 
tested the diffraction caused by a single vortex imposed 
on the background beam. As seen in Fig. 6e the diffrac- 
tion by a single vortex is substantially different and re- 
sembles the diffraction of a laser beam by a wire. In 
all our experiments we observed strong vertical asymme- 
try of the probe-beam diffraction pattern, which always 
developed downwards at higher powers. Numerical mod- 
eling of the processes and further experimental investiga- 
tions should allow to gain deeper insight in the relative 
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strength of the mentioned mechanisms. 



V. CONCLUSION 

In conclusion, we successfully generated experimen- 
tally lattice structures of optical vortices with different 
topological charge distribution and described their prop- 
agation in saturable NLM. Due to the intensity depen- 
dence of the refractive index these lattices induce periodic 
modulation of the refractive index of the medium and 
'write' and effective phase grating in it. The modulation 
is sufficient to force a perpendicularly propagating probe 
beam of a He-Ne laser to diffract. This feature could ap- 
pear as an interesting possibility to create periodic struc- 
tures in the refractive index of a NLM. It could find an 
application for optical writing of two-dimensional pho- 
tonic crystals and could appear relevant to the physics of 
Bose-Einstein condensates. 
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